The Virial Theorem in Graphene and other Dirac Materials 
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The virial theorem is applied to graphene and other Dirac Materials for systems close to the 
Dirac points where the dispersion relation is linear. From this, we find the exact form for the total 
energy given by E — B/r s where r s oo is the mean radius of the d-dimensional sphere containing one 
particle, with ao the Bohr radius, and B is a constant independent of r s . This result implies that, 
given a linear dispersion and a Coulombic interaction, there is no Wigner crystalization and that 
calculating B or measuring at any value of r s determines the energy and compressibility for all r 3 . 
In addition to the total energy we calculate the exact forms of the chemical potential, pressure and 
inverse compressibility in arbitrary dimension. 
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The linear dispersion of graphene near the Dirac 
points 1] combined with the experimental realization of 
graphene sheets [2| has provided the impetus for count- 
less theoretical and experimental investigations [3] . Stud- 
ies of graphene have been fueled by the material's broad 
spectrum of interesting features — from the pragmatic 
possible technological implementations to its relationship 
to massless fermions travelling at the speed of light, as the 
linear dispersion is identical in form to the Dirac Hamil- 
tonian for massless relativistic fermions. From the point 
of fundamental physics, one of the most important ques- 
tions is about the role of electron-electron interaction in 
the system with linear dispersion. This question is more 
complicated than it sounds for the Dirac Materials, ma- 
terials that have linearly dispersive excitations be- 
cause at least close to the Dirac point, where the Fermi- 
surface reduces to a point, the density of state vanishes 
giving rise to a suppressed Coulomb interaction screen- 
ing. For reduced density, the Coulomb interaction in con- 
ventional electron systems would lead to the localization 
of the charge in space, also known as Wigner crystal- 
lization. Does this conventional physics hold in linearly 
dispersive system at the dilute limit? or does the differ- 
ence in the dispersion relation introduces unconventional 
physics? This question is very relevant because the local- 
ization of carriers in graphene has not yet been observed. 

Central to the determination of the role of the Coulomb 
interaction is the measurement and calculation of the 
electronic compressibility which depends upon the total 
of the potential and kinetic energy. The compressibil- 
ity of monolayer graphene has been measured and the 
data seems to be well described by the kinetic energy's 
contribution alone || with effective Fermi velocity which 
incorporates the effect of the Coulomb interaction. Ran- 
dom Phase Approximation [f| and Hartree-Fock [?J stud- 



ies have concluded that contributions to the compress- 
ibility from the Coulomb interaction yield a (logarithmi- 
cally) doping dependent correction of between 10 and 20 
percent in the experimental regime. We see significantly 
different phenomenon viewing the same problem through 
the lens of the virial theorem. 

Using the virial theorem for the linearly dispersive sys- 
tem, we find that an exact form of the total energy valid 
for all doping densities can be expressed as 



B 

E{r s ) = f 



(1) 



B is a constant which does not depend on r s , however, 
the exact numerical value will depend on the material 
and the dimensionality of the system. Eq. (flj is the 
main result of this work. The virial also gives concise 
confirmation of the absence of Wigner crystalization at 
all densities originally argued in Q for low densities, and 
a qualitatively different description of the compressibility 
than is currently reported in the literature. 

The approach we use to obtain Eq. |T]) follows the argu- 
ment of Argyres [9[ , applied earlier to a quadratic disper- 
sion, with slight but substantive modifications to account 
for the linear dispersion. We consider the virial theorem 
for a system with a single particle energy that has a linear 
dispersion relation E p = ±v • p, with v = v g p (where in 
graphene the characteristic velocity v g ~ 10 6 m/s), and 
interacts with the other electrons and a uniform back- 
ground charge via a Coulombic potential V oc 1/r. It is 
in this sense that this result is valid for other (Coulombic) 
Dirac materials found in arbitrary dimension. The virial 
has been studied extensively for the conventional electron 
gas 0, but has not yet been studied in graphene or 
other Dirac materials with a linear dispersion. 

We focus here on the two dimensional case and intro- 
duce the generalized volume £1, which is the area A in 
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2D, but the proof applies in arbitrary dimension. Given 
a Hamiltonian H = T + V , where the kinetic energy 
is given by T — • Pi, and the potential energy 

V = V e + V e b + V;,. V e is the potential energy between 
electrons, 



V K = 



-y- 



(2) 



V e b is the potential energy between the electrons and the 
background of positive charge, 



d 2 x- 



(3) 



and Vb is the potential energy of the background of pos- 
itive charge, 



V b = I ( d 2 x 

Z J A J A 



d 2 x'^- 



(4) 



The background positive charge density p oc r~ 2 and the 
area Aar,. V e is then a homogeneous function of degree 
-1 in {r^}, Vb is a homogeneous function of degree -1 in 
r s , and V e b is a homogeneous function of degree -1 in 
{ri,r s }. The Euler theorem for homogeneous functions 
combined with the definition of V yields 
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dV dV 
Yi ■ = — V — r s . 
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(5) 



Eq.[5]is exactly the same result one gets when calculating 
the virial for the 2DEG. 

The difference lies in the kinetic term, for the linear 
(quadratic) dispersion in graphene (2DEG) the kinetic 
energy is a homogeneous function of degree 1 (2). Thus, 
for graphene 



x dT 

i d & 
We consider the virial 

G = r i ■ P« 
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the quantum mechanical time derivative of which is 

G = hG,H] 
in 



(G) 



(7) 
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The last equality follows from Eqs. §5§ and (J6)). Because 
T is independent of r s , r. 



dv 



dH 



expectation value of Eq. ([8]) , note that [G, H] 
find El 



We then take the 
and 



T + V 
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(9) 



The remarkable feature of Eq. which is not present 
in its 2DEG equivalent with parabolic dispersion is that 
the left hand side is simply the average total energy, such 
that 

E = -r a ^. (10) 

Thus, the average total energy for many-body systems 
with a linear dispersion interacting via a Coulombic po- 
tential satisfies a first order differential equation; the so- 
lution, Eq. ([1]), has many implications including those 
discussed below [26j . Eq. (fTU| is valid for all d, but other 
thermodynamic properties will have some explicit depen- 
dence on the dimension. 



Following the work of Ceperley, as in , it can easily 
be shown that 

where P is the pressure, fl is the generalized volume and 
d is the dimension. For a linear dispersion, we have 

T+V _ E 



pn = 



(12) 



where again we have a simplification leading to a direct 
relationship involving the total energy E and, in this case, 
thermodynamic quantities. This can be readilly used in 
a calculation of the chemical potential for graphene as 
well as other Dirac materials in arbitrary dimension via 
the Gibbs free energy 

1 + 1 ) E = fiN, (13) 



G = E + Pfl 



so. 
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(14) 



From the expression for the total energy, Eq. (1), we 
can derive the electronic compressibility k for T = 0, 
constant N and d = 1, 2 and 3. We will focus first on the 
compressibility for 2D since this case, with a spin degen- 
eracy of 2 and valley degeneracy of 2, models the Dirac 
material graphene. To compare the theoretical predic- 
tion of k with the measurements done on graphene [2], 
it is convenient to work with the energy per unit area, 
£ = E/A, we then find that 

8 2 £ 
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(15) 
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where we have introduced the expression k — n 2 K and 
n = N/A. Our result, Eq. (|15|) . which follows directly 



from the Virial Theorem, shows that R 



-1/2 



cx n^ 1 -'", as 

seen in experiment [5|, since B is independent of n. In 
their paper, Martin, et al. [5| have shown that they can 
fit their data using the 2D free Fermi gas with a Dirac 
spectrum and a spin + valley degeneracy of 4. In this 
limit their only fitting parameter is the velocity, v g , where 
v g ~ 10 6 m/s is close to the value obtained from the band 
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structure 

While this interpretation is consistent with Eq. (fT"5j) we 
need to look more closely at the interpretation of Eq. (fl"5|) 
from the point of view of the 2D Fermi liquid theory to 
get a better understanding of the role of interactions in 
the Dirac materials. The derivation of h in a 2D Fermi 
liquid with a Dirac spectrum is a straightforward gener- 
alization of the 3D Landau Fermi liquid result [12J, LL3J , 
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(16) 



with 



where the density of states is given by iV(0) = 

Vp the renormalized Fermi velocity and Fq is the spin 
symmetric quasiparticle interaction 12J. The Fermi mo- 
mentum is related to the density by n = kp/ir, thus 
(1 + Fq )v* f = (3B/2N)a must be independent of the 
density to be consistent with Eq. (fPo)) . When compar- 
ing Eqs. (15) and (16) it is clear that only measuring k 
is not sufficient to determine the renormalizations of v F 
and 1 + Fq independently. To determine the strength 
of the interaction parameter Fq we need an independent 
measurement of the density of states. This can be ob- 
tained from the low temperature electronic specific heat 
at constant volume, 



c v = —N(0)T. 



(17) 



With a measurement of the specific heat we can an- 
swer the question: Do many body effects, coming from 
electron-electron interactions play an important role in 
graphene? 

In lieu of the measurements, we can make some com- 
ments regarding interaction effects in graphene. In refer- 
ence Q Dahal et al. argued that Wigner crystallization 
could not occur in graphene based on the Hartree-Fock 
results for the kinetic T/N and potential V/N energy 
per particle, where 
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(18) 



e is the charge of the electron and e is the dielectric 
of the medium calculated in the RPA [S] approximation 
and is independent of r s . They found [8j that the ratio 
T/|y| ~ 1 is independent of r s . Thus, the coefficient B 
calculated in the H-F approximation, B HF 



is positive and independent of r s but is much smaller than 
the experimental value To improve on the H-F ap- 
proximation for S HF , one needs to include higher order 
corrections to the ground state energy, for example, the 
RPA corrections used in reference |6J . Alternatively, one 
could do the H-F approximation for the self-energy 0], 
which gives rise to a renormalization of the Fermi velocity 
v F . While the approaches used in @ and Q, are clearly 
the right thing to do, the end result is, both approaches 
give rise to a ln(l/r s ) contribution to the compressibility 
that is not seen experimentally. 

The approximations, H-F and RPA, used in refer- 
ences 0, [3] and [1] for calculating the energy and the 
compressibility have worked well as starting points in the 
traditional 2DEG and 3DEG. The reasons for this are: 
1) The H-F approximation is considered to be exact in 
the low-density limit (r s 3> 1) and 2) The RPA is con- 
sidered to be exact for the high-density limit (r s <C 1). 
Thus, in the traditional Coulomb gas, we have regimes 
of density where we can add perturbative corrections to 
these leading order exact results. From reference and 
more precisely from Eq. ([!} it is clear that there is no 
perturbative regime for graphene, or more generally for 
the case of a charged 2 dimensional Dirac gas (2DDG), 
since T/lV"! ~ 1 for all r s . What we have learned from 
the virial theorem is that it is not surprising that the 
perturbation theory based efforts to calculate the 

compressibility have not agreed with experiment. Most 
likely calculating the energy and/or the compressibility 
in a 2DDG would require a numerical calculation or a 
non-perturbative many-body approach, along the lines 
ofRef. Q. 

The effects of the interactions are even more striking 
when we consider the 1 dimensional Dirac gas (1DDG). 
This is not just a theoretical toy model since it is pos- 
sible to make a 1DDG by rolling up a graphene sheet 
along specific directions, forming a metallic, ID, single 
walled carbon nanotube [r|. There is a lot that is 
known about interac ting Fermions in ID, see for example 
reviews in Refs. 17j,|2l|. One of the more important re- 
sults is that in ID the interactions between the fermions, 
no matter how weak, turn a Fermi gas into a Luttingcr 
liquid, see for example Refs. (l8l - [2l| . Another important 
result is that in the case of the 1DEG there is no Wigner 
crystal formation at any density [22J. This is expected 
since a long range ordered state, like a Wigner crystal, 
cannot exist in ID [2jJ. 

The absence of the Wigner crystal can be seen more 
clearly in terms of the energy density £ = E/L, where E 
is given by Eq. ([I]), 
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We have defined n = N/L, with L the length of the ID 
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chain, so n = and n — Akp/n. Eq. (| 19b|) for k is 
most interesting since it is independent of the density for 
ID. We expect this behavior to be seen in single walled 
carbon nanotubes. It turns out that there are models in 
ID that satisfy Eq. (|19bp , and one of them, we will see, is 
the Tomonaga-Luttinger model. In this model there is an 
exact solution for the limit of extremely long-range inter- 
actions when the non-interacting single particle spectrum 
is lincraized. This problem can be solved by the standard 
bosonization techniques, see, for example Ref. (2ll |. turn- 
ing the interacting fermions into non-interacting bosons 

The boson velocity, i.e., the speed of 1st sound, ci, is 
given by, see [13], c 2 = v f + (2vf/tt)\, where vf — kp/m 
and to is the bare electron mass. Here the electron- 
electron interaction in momentum space is approximated 
by a coupling constant A which is taken to be indepen- 
dent of the density and drops off rapidly for momentum 
transfers q <C kp. If we plug this expression for c\ [2(| 
into the compressibility k^ 1 = (m/n)c{, we get a k~ l 
that depends on the density. If we instead use the Dirac 
spectrum where vf is replaced by v g and kp/v g replaces 
to, we get k _1 = (fcj? /v g n)c\. We find a density indepen- 
dent compressibility and first sound, c 2 ; we also find an 
energy density that is proportional to n 2 [28| . 

From the Tomonaga-Luttinger model, there is an in- 
sight here that is relevant to the 2DDG. If we calculate 
the energy density and compressibility for the ID free 
Dirac gas we would get that £ cx v g n 2 and k" 1 cx v g , 
which has the same density dependence as the interact- 
ing model. Clearly, we cannot conclude that the inter- 
actions are not important simply because we can get the 
same result as the Tomonaga-Luttinger model by simply 
rescaling the Fermi velocity in the non-interacting Dirac 
gas. The point here is that in the 1DDG, no matter how 
weak the interactions, the ground state of this model is 
not a renormalized Fermi liquid, but in fact a Luttinger 
liquid. 

Furter evidence of the power of the virial theorem in 
Dirac materials can be seen completely outside of con- 
densed matter, via studying weakly interacting quarks in 
the massless limit as in the MIT bag model [23]. Baym 
and Chin [24( find that to order g 2 in the color coupling 
constant g c , the speed of sound c\ is independent of den- 
sity, precisely as predicted by the virial theorem for a 
linear dispersion in 3D. 

In this letter we have presented the virial theorem for 
linear kinetic and Coulombic potential energies, yielding 
a simple differential equation for the total energy, and we 
predict that E ~ l/r s for all densities in Dirac Materials. 
We have presented the wide applicability of this power- 
ful new solution in 1, 2 and 3 dimensions and it has been 
shown to agree with the experiment for the compressibil- 
ity in 2D over a wide density range. We have suggested 
that the results required by the virial theorem for the 
linear spectrum with a 1/r potential apply to other po- 



tentials and presented systems in which this is the case. 
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